Is the Dynamics of Scaling Dark Energy Detectable? 
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We highlight the unexpected impact of nucleosynthesis and other early universe constraints on 
the detectability of scaling quintessence dynamics at late times, showing that such dynamics may 
well be invisible until the unveiling of the Stage- IV dark energy experiments (DUNE, JDEM, LSST, 
SKA). Nucleosynthesis strongly limits potential deviations from ACDM. Surprisingly, the standard 
Chevallier-Polarski-Linder (CPL) parametrisation, w(z) = wq + w a z/(l + z), cannot match the 
nucleosynthesis bound for minimally coupled scaling fields. Given that such models are arguably 
the best-motivated alternatives to a cosmological constant these results may significantly impact 
future cosmological survey design and imply that dark energy may well be dynamical even if we do 
not detect any dynamics in the next decade. 



I. INTRODUCTION 

The quintessential enigma of modern cosmology is that 
of dark energy - the shadowy source of the accelerated 
expansion of the cosmos. Einstein's "greatest blunder" - 
the cosmological constant, A, remains the most econom- 
ical explanation for this observed acceleration, but more 
exciting alternatives could imply the existence of new 
types of matter, modifications of the Einstein equations 
or even violation of the Copernican principle. 

Of these however, one can argue that only scalar fields 
are consistent with current constraints while simultane- 
ously being theoretically well-founded. The more ex- 
otic alternatives either are claimed to have theoretical 
pathologies or are likely to be indistinguishable from the 
vanilla ACDM model which, given all current data, re- 
mains flavour of the month [l(. Models in the former 
category include K-essence Q which may exhibit an un- 
bounded adiabatic speed of sound Q and the DGP model 
[H, which has some theoretical problems with ghosts 
(see Q and references therein) while those in the lat- 
ter class include f(R) modifications of gravity @ and 
unified dark energy [3]. Of the scalar fields, arguably 
the best-motivated and most compelling are the scaling 
quintessence models [1, H [l(| ■ We consider such scal- 
ing models with a redshift-dependent equation of state 
parameter w(z) — p/p that tracks the dominant energy 
density component of the cosmos (w = | , 0) until a red- 
shift z = Zt at which point it undergoes a transition to 
w < which triggers the onset of acceleration. 

Big Bang Nucleosynthesis (BBN) provides strong con- 
straints on the energy density of dark energy during the 
radiation dominated era at a temperature of T ~ lMeV, 
implying that n DE (T ~ lMeV) < e = 0.045 at 2a [HI- If 
the scalar field has reached the scaling attractor solution 
at the time of BBN, this early universe constraint on the 
energy density of the dark energy will be carried forward 
until the field exists the scaling regime (at z t ). The time 



at which the attractor for the scalar field is approached 
depends on the details of the physics before the radiation 
dominated epoch. But in [ljj] it was shown that, for a 
typical inflationary model with the usual method of re- 
heating, the field will approach the attractor long before 
nucleosynthesis. Therefore we assume that the field is 
already scaling with the background fluid at the time of 
BBN. In fact, our results are unchanged even if scaling 
only occurs by decoupling, since similar or better con- 
straints on the dark energy density exist from the Cosmic 
Microwave Background (CMB), fi DE (T ~ leV) < 0.04 
12]. Constraints on w(z) from BBN taken together with 
recent data require that z t > 5 [13| . 

In the case of scaling dark energy models these con- 
straints at early times are preserved by the subsequent 
evolution of the cosmos, implying that the limit Qde < 
e = 0.045 holds for all redshifts greater than z t , the red- 
shift at which the equation of state stops scaling and 
becomes negative. We find new implications for the mag- 
nitude of allowed deviations of these scaling models from 
A. To motivate our results consider first a toy, step- 
function model for w(z) in which w is constant for z < z t 
and w = for z > z t . Requiring OdeO^) < e = 0.045 
(the BBN constraint preserved by the scaling field) with 
w > — 1 and f^DE = 0.7 today implies that Zt > 2.6, 
an unexpectedly large number given that most cosmic 
probes in the next decade will be limited to z < 2. 
To achieve z t — 1 one requires instead w = —1.88, a 
value disfavoured by current observations [l[ . Figure © 
shows that the early dark energy constraints on the scal- 
ing models we consider force the derivative of w(z) to be 
very small: |u/(0)| < 0.2 for all our models, significantly 
smaller than will be detectable in the next decade [15j . If 
w(0) < —0.9 the deviation of the Hubble rate from A is 
less than 3% and the deviation of the distance modulus 
is less than 0.04 mag. In addition, we find the surprising 
result that the standard CPL parametrisation [lj] can- 
not be used to describe scaling minimally coupled scalar 
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z t ) — 0, one can show that 

. e(l + z) 3 0.047, 

f(z > z t ) = = 1 + z 3 

r(l — e) r 



(2) 



where r — Q,w/Q, m is evaluated today. 

Using this result implies the following general but 
stringent constraint: 
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= 1.023. 
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FIG. 1: Evolution of 17de(z) for the models we consider, 
showing their approach to the BBN limits of e — Qve{z = 
zbbn) = 0.045 for the polynomial parameterisation of w{z) 
and |e for the double exponential potential. For comparison 
we also show the curves for A and the CPL [3] w(z) with 
the lowest asymptotic value of S7de in this model while still 
assuming w > — 1, showing its inability to match the BBN 
constraint and describe a canonical scalar field. Figure ((2]) 
shows the corresponding observational quantities for the scal- 
ing quintessence models. 



fields (which automatically have w > — 1) and simultane- 
ously match the nucleosynthesis bound. Throughout the 
paper we assume a flat universe with f2 m = 0.3 today 
and all results are valid for redshifts less than matter- 
radiation equality. We start by describing the constraints 
on general models in Section [Til we gi ye the specific scal- 
ing models we consider in Section IIIII and discuss the 
performance of standard parameterisations in describing 
these scaling fields in Section [iVl We conclude in Section 

El 



II. GENERAL RESULTS 

One can derive model-independent constraints on cos- 
mological observables (we consider the Hubble rate, 
H(z), and distance modulus, (J,(z)) for scaling models. 
The BBN constraint implies that ^de(z > z%) < e = 
0.045 since J7de(z > z t ) is constant in perfectly scaling 
models. It is standard to write pbe{z) = HqQ,^(z = 
0)f(z) where 



f{z) = exp 



1 + w(z') 
1 + z' 



dz' 



(1) 



determines the redshift dependence of the dark energy 
density. The evolution of Ode as a function of redshift 
is illustrated for various models in Figure ([T]). From the 
above discussion, plus the scaling requirement that w(z > 



where the last equality arises from imposing e = 0.045 
and the upper bound from setting r = as z t — > oo. 
This limit can clearly be seen in the top middle panel 
of Fig. ([2]). This robust result implies that detecting 
deviations from A at high redshift will be difficult using 
Hubble rate measurements alone since they are bound 
to be less than 2.3%. This does not, however, strongly 
constrain the behaviour of Hde/Ha for z < z t , once the 
field has left the scaling regime. We will show in two 
classes of models that its maximum value is less than 5% 
and occurs around z ~ 1. The latter result is good news 
for Baryon Acoustic Oscillation (BAO) surveys such as 
WiggleZ, BOSS and WFMOS [13] which will probe this 
range of redshifts. 

Similarly one can also place robust bounds on the de- 
viation of the distance modulus, fJ,(z), from the A predic- 
tion. The quantity A/j, = ^de(z) — Ma(z) is given by: 



A^i = 5 log 



10 



4,de(^) 



(4) 



where cIl(z) is the corresponding model luminosity dis- 
tance. If we assume that there exists a number a such 
that H(z)j H\(z) < 1 + a 2 for all z (for example in 
the constraint on H discussed above a 2 ~ 2.3%) then 
dL,DE{z)/dL : \{z) > (1 + a 2 ) -1 and hence A/i(z) obeys 
the inequality 



0> Afi(z) > -51og 10 (l + a 2 ) 



(5) 



A bound of a 2 = 0.025 gives |A/i(z)| < 0.054 mag - see 
the top row of Fig. This is a conservative upper 

bound since we have shown that for z > Zt, ot 2 < e/2. 
For the oscillating double exponential w(z) models (with 
—30 < (i < 0) considered later, one has a 2 < 0.015 which 
yields the constraint |A/z(z)| < 0.032 mag. These general 
results do not constrain H(z) for z < z t , which instead 
requires a specific model for w(z). We now consider two 
classes of models which describe a wide range of scalar 
field dynamics. 
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FIG. 2: BBN compatible scaling dark energy models: allowed observables for the polynomial w(z) for z t from 4.02 (light 
brown) to 10.02 (dark brown) (top panel) and double exponential potential with —30 < fi < 1 (bottom) showing the small 
deviations from the A predictions. The bottom light brown curves are from fi = 1 (for which w(0) = —0.8) to fi = 0, w(0) = — 1 
and the dark brown curves are for fi = -30, -25, -20, -15, -10, -5, -4, -3, -2, -1, -0.8, -0.6, -0.4, -0.2. All models have 
\dw/dz(0)\ < 0.2. In the model with w(0) < —0.9 Hue(z) deviates from Ha{z) by at most 2.7% (marked by the horizontal line 
in the bottom middle panel) and implies a deviation in distance modulus of less than 0.05 mag. The error bars in the right 
panels correspond to the Stage-Ill (large boxed errors) and Stage-IV (small triangular errors) supernova surveys respectively. 
These are produced for the bottom curve in each case except for the Stage-IV (SNAP-like) errors in the double exponential 
case which correspond to the w(0) = —0.9 model (thicker line asymptoting to ~ 0.04 mag). Note that the ratio H(z)/Ha for 
the double exponential potential does not converge to 1.023 since the matter-dominated value of Qt>e{z) is constrained to be 
3/4 of the radiation-dominated value. Hence the Hubble rate is forced to ~ 1 + (3/8)e ~ 1.017. 



III. SPECIFIC FORMS OF THE SCALING 
POTENTIAL 

A. Polynomial w(z) parametrisation 

First we consider a quadratic parametrisation of the 
dark energy equation of state, w(z) [17] : 

i \ I w o + w i z + W2Z 2 for z < Zt ,„« 



We apply the constraint w{z) > —1 since we want to 
describe minimally coupled scalar fields with canonical 
kinetic terms. The linear case with W2 = requires z t ~ 
6.2 to match BBN (for e = 0.045) if w a = -1 and if 
we allow wq to be free the BBN constraint implies the 
correlation w\ ~ —OAwq — 0.2 for the interesting region 
— 1 < wq < —0.8. The other case we consider is wq = 



—1,W2 =/= 0. Continuity at z — z t then implies 

U >2 = -2-—- ( 7 ) 

z t Z t 

The BBN constraint provides w\ in terms of z t . The re- 
sulting family of curves and observables are shown in the 
top row of Fig. ([2]). We note that for z < 1 the BBN con- 
straint is so strong that all models have w(z) < —0.8, as 
shown in Fig. ((2]), and the largest deviation of H(z) from 
the A comparison model is about 2.7% with the largest 
deviation in the distance modulus only about 0.03 magni- 
tudes (occurring at z = 2). This shows that if w ~ —1 to- 
day we cannot expect significant deviations from A at any 
redshift and only Stage-IV experiments [lj| are likely to 
detect dark energy dynamics with any real significance. 
However while current data favour w(0) ~ —1 they are 
consistent with larger values. To study how this zero- 
point affects our results we now consider simulations of 
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a scalar field with a popular family of scalar field poten- 
tials, V (<f>), which also allows us to examine the impact 
of oscillations in w{z). 



B. Double Exponential Potential 

A single exponential potential is well known to give 
early scaling Q but cannot also lead to late-time accel- 
eration. One well-studied way to combine the two is via 
the double exponential potential [Tj|: 



V{4>) = Mfe-^ + M 2 V 



flK-4 



(8) 



where late-time acceleration is induced if p < The 
dynamics of (f> is sensitive to the sign of p, in that nega- 
tive values of p yield a potential that oscillates around the 
minimum of the field. This can clearly be seen in Fig. [21 
for the potentials with — 30 < p < We assume that 
the scalar field is in the scaling regime during radiation 
domination, and hence satisfying the BBN constraint re- 
quires A > 2/Ve > 9.43 [g. We choose A = 9.43 to 
maximise deviations from A. Unlike a perfectly scaling 
model of the form we assumed in the previous section, 
Qcf, actually decreases in the transition to matter domi- 
nation and we have fl^ < 3/4 x e = 0.034 during matter 
domination, leading to even more stringent results - see 
Fig. ([I]). The extra 3/4 factor is specific to the double 
exponential potential and is responsible for the reduc- 
tion of the asymptotic values of Hw(z) / Hacdm, seen in 
the bottom middle panel of Fig. ((2]), relative to the pre- 
dictions of Eq.([3]). We choose Mi = 10~ 14 m pl and use 
Planck units where k = 1. 

We numerically solve the evolution equation 



3H <f> + Va, = 



(9) 



for <fi with V(4>) given by Eq.©, together with the radi- 
ation and matter fluids. These are all coupled to gravity 
through the Fricdmann equation 



H 2 = — 
3 



-ft + V{4>)+ Pr + Pn 



(10) 



For each value of p we find M-i such that J7de = 0.7 today 
which implies M 2 ~ 10~ 31 m p i for p ~ 1. The resulting 
w{z) curves for these models, together with the predicted 
observables (H(z), Ap(z)) are shown in the bottom row 
of Fig. @. For negative p - we study —30 < p < - 
the potential has a global minimum and for z < 0.2 the 
equation of state satisfies < —0.98. As a result all 
the negative p models show tiny deviations from A: less 
than 1.5% for H(z) and less than 0.015 mag for Ap (see 
Fig. $2$). This will make detection extremely difficult 
even with the Stage-IV dark energy experiments such as 
DUNE, JDEM, LSST and SKA [lj|. 

In contrast, positive values of p (0 < p < 1) simply 
modify the slope of the potential. They can yield values 
of significantly different from —1 today, e.g. for p ~ 1 



one finds w^0) ~ —0.8 which is consistent (at about the 
2cr level P, [l3j]) with current observations and which we 
therefore take as the upper bound for p. The bottom 
row of Fig. @ shows that the maximum allowed devia- 
tion for H(z) from H/^(z) in this case is about 5%, peak- 
ing at z ~ 1 with a maximum value of Ap ~ 0.9 mag. 
Such a model will be detectable with Stage-Ill supernova 
surveys (at the 99.97% confidence level) and with the up- 
coming BAO experiments, since the maximum deviation 
in H(z) coincides with the redshift ranges in which they 
will operate, i.e. z ~ 0.7 to 1.1. However, for values 
more consistent with the current best-fits, w^O) < —0.9 
one finds much smaller deviations of 2.7% and 0.045 mag 
respectively for H(z) and Ap which again will require 
Stage-IV experiments for conclusive detection as can be 
seen in the right-hand panels of Fig. @. Similar results 
will apply to other modifications of the exponential po- 
tential, e.g. fl9l ]. 



IV. PERFORMANCE OF STANDARD 
PARAMETRISATIONS 

The most widely used parametrisation for dark en- 
ergy, the Chcvallicr-Polarski-Linder (CPL) parametrisa- 
tion G3, 



w{z) = Wq + 1 



l + z' 



(11) 



is also the basis for the DETF figure of merit [15j. Sur- 
prisingly this parametrisation fails dramatically to meet 
the BBN constraint if we demand w(z) > —1, £!de ~ 0.7 
today and w(z > z t ) — as before (for some z t ). This is 
clearly visible in the top curve of Fig. ((T|), which shows 
the lowest attainable value of SIde with wq = — 1. In ret- 
rospect this is understandable since to reach the scaling 
value w = for some z t requires w a > —wq. However in 
this case w(z) doesn't spend enough time at sufficiently 
negative values to force JIde^) down to the BBN value. 
The least phantom value of Wq that is able to satisfy the 
BBN constraint is wo = —1.3. In contrast the logarithmic 
expansion w(z) = wq + w± ln(l + z) is able to match the 
BBN constraint with w{z) > —1, but only for z t > 12.4. 



V. CONCLUSIONS 

Scaling field models are arguably the best-motivated 
alternatives to the cosmological constant. We have shown 
that the constraints on the energy density of the scalar 
field at the time of Big Bang Nucleosynthesis and decou- 
pling strongly limit their allowed dynamics. If w today 
is not close to the maximum value allowed by current 
data then detection of dynamics will likely have to wait 
a decade for the Stage-IV DETF experiments. Of course, 
these strong conclusions are only true for scaling models 
and if one allows exotic phantom behaviour (w < — 1) 
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the conclusion is much more rosy. We discussed two spe- 
cific families of these scaling models - imposing the BBN 
and CMB constraints on general scaling models is left for 
future work. 

One might ask how the constraints discussed here will 
affect growth of structure. Since the growth in minimally 
coupled scalar fields smoothly approaches that of ACDM 
when w(z) — > — 1 (see e.g. Fig. 2 in [20() our results 
show that the growth in models allowed by BBN and 
CMB constraints will be close to those in ACDM at z < 
1, making detection difficult. To what extent there are 
deviations at higher redshift, and whether they will be 
detectable by future cluster or weak lensing surveys is 
again left as a subject for future study. 

Finally we have shown that the standard CPL 
parametrisation, w(z) = wq + w a z/(l + z), fails dra- 
matically to match the BBN constraint when describing 
scaling fields which satisfy w > — 1. This is particularly 
important given that the CPL parametrisation is the ba- 



sis of the DETF figure of merit [lj| which is now the 
de facto standard for the optimisation of future cosmo- 
logical surveys, e.g. [2l[. A concern therefore is that 
optimisations may be unwittingly biased away from scal- 
ing dark energy models. More work in this area is clearly 
needed to assess the implications for cosmological survey 
design, but it is clear that the current non-detection of 
dark energy dynamics should neither come as a surprise, 
nor should it discourage us from the hunt. 
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